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■ Abstract 

in 

Orthogonal Frequency Division Multiplexing (OFDM) has been combined with Multi-Input Multi-Output 
(MIMO) techniques for broadband wireless communication systems. Bit-Interleaved Coded Multiple Beamforming 
\ (BICMB) can achieve both spatial diversity and spatial multiplexing for MIMO flat fading channels. For MIMO 

frequency selective channels, BICMB with OFDM (BICMB-OFDM) can be applied to achieve both spatial diversity 
and multipath diversity, making it an important technique. However, analyzing the diversity of BICMB-OFDM is 

> 

a challenging problem. In this paper, the diversity analysis of BICMB-OFDM is carried out. First, the maximum 
" achievable diversity is derived and a full diversity condition R C SL < 1 is proved, where R c , S, and L are the code 

en; 

. rate, the number of parallel streams transmitted at each subcarrier, and the number of channel taps, respectively. 

Then, the performance degradation due to the correlation among subcarriers is investigated. Finally, the subcarrier 
grouping technique is employed to combat the performance degradation and provide multi-user compatibility. 
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I. INTRODUCTION 

Substantial research and development interest have been drawn on Multiple-Input Multiple-Output 
(MIMO) systems because they provide high spectral efficiency and performance in a given bandwidth. In 
a MIMO system, beamforming techniques exploiting Singular Value Decomposition (SVD) are employed 
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to achieve spatial multiplexing^ and thereby increase the data rate, or to enhance performance, when the 
Channel State Information (CSI) is available at both the transmitter and receiver 

For flat fading channels, single beamforming carrying only one symbol at a time achieves full diversity 
01, 0. However, spatial multiplexing without channel coding results in the loss of the full diversity order. 
To overcome the performance degradation, Bit-Interleaved Coded Multiple Beamforming (BICMB) was 
proposed J6]|, 0, [H. BICMB systems studied so far employ convolutional codes [|9l as channel coding, 
and interleave the coded bit codewords through the multiple subchannels with different diversity orders. 
BICMB can achieve the full diversity order as long as the code rate R c and the number of employed 
subchannels S satisfy the condition R C S < 1 IfTOl , IfTTl . 

If the channel is frequency selective, Orthogonal Frequency Division Multiplexing (OFDM) can be used 
to combat the Inter-Symbol Interference (ISI) caused by multipath propagation [fT2l . OFDM employs a 
large number of closely spaced orthogonal subcarrier signals to carry data, which are divided into several 
parallel data streams. In addition to its robustness against ISI, OFDM achieves many advantages such as 
easy adaptation to severe channel conditions without the computational burden of time-domain equaliza- 
tion, high spectral efficiency, efficient implementation using Fast Fourier Transform (FFT) and Inverse FFT 
(IFFT), and low sensitivity to time synchronization errors. Furthermore, multipath diversity can be achieved 
by adding channel coding |[T3l . |fl4l . OFDM is well-suited for broadband data transmission, and it has 
been selected as the air interface for the Institute of Electrical and Electronics Engineers (IEEE) 802.11 
Wireless Fidelity (WiFi) standard, the IEEE 802.16 Worldwide Interoperability for Microwave Access 
(WiMAX) standard, as well as the Third Generation Partnership Project (3GPP) Long Term Evolution 
(LTE) standard lfT5H . 

MIMO techniques have been incorporated with OFDM for all broadband wireless communication 
standards, i.e., WiFi IfToll , WiMax [fTTl , and LTE [fT8l . Beamforming can be combined with OFDM for 
frequency selective MIMO channels to combat ISI and achieve spatial diversity |fT9l . Moreover, both 
spatial diversity and multipath diversity can be achieved by adding channel coding, e.g., BICMB with 
OFDM (BICMB-OFDM), [|20]|, Efl, 0. Although some of the more modern standards employ more 
sophisticated codes lfT5l . such as turbo codes and Low-Density Parity-Check (LDPC) codes [9], than 
convolutional codes employed by BICMB-OFDM, convolutional codes are important since they can be 
analyzed and there is a great deal of legacy products using them. Therefore, BICMB-OFDM is an important 

'in this paper, the term "spatial multiplexing" is used to describe the number of spatial subchannels, as in (T|. Note that the term is 
different from "spatial multiplexing gain" defined in (2)- 
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Fig. 1. Structure of BICMB. 



technique for broadband wireless communication. However, the diversity analysis of BICMB-OFDM is a 
difficult challenge. 

In this paper, the diversity analysis of BICMB-OFDM is carried out. First, the maximum achievable 
diversity is derived and the important a-spectra directly determining the diversity are introduced. Based 
on the analysis, a sufficient and necessary full diversity condition, R C SL < 1 where S is the number 
of steams transmitted at each subcarrier and L is the number of channel taps, is proved. Then, the 
performance degradation caused by the correlation among subcarriers is investigated. To overcome the 
performance degradation, the subcarrier grouping technique ll22l . Il23l . Il24ll which also provides multi-user 
compatibility, is employed. 

The remainder of this paper is organized as follows: Section |n] briefly describes BICMB and OFDM. 
In Section Unl the system model of BICMB-OFDM is outlined. In Section [TV] the maximum achievable 
diversity is derived, an important a-spectrum directly determining the diversity is introduced, and a 
sufficient and necessary condition to achieve full diversity is proved. Then, Section |V] investigates the 
performance degradation caused by the correlation among subcarriers and Section [VI] discusses the 
subcarrier grouping technique to combat the performance degradation and provide multi-user compatibility. 
In Section IVIIl simulation results are provided. Finally, a conclusion is provided in Section IVIIII 

II. BICMB AND OFDM 

A. BICMB 

For flat fading MIMO channels, the diversity of SVD beamforming is 

D = (N r -S + l)(N t -S + l), (1) 

where N r and N t denote the number of transmit and receive antennas respectively and S < min{N t , N r } 
denote the number of streams transmitted at the same time [4J. BICMB was proposed to overcome the 
loss of full diversity N r N t for spatial multiplexing, i.e, S > 1 0, 0. i- 
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The structure of BICMB is presented in Fig. [TJ First, the convolutional encoder of code rate R c , possibly 
combined with a perforation matrix for a high rate punctured code Il25l , generates the bit codeword c 
from the information bits b. Then, an interleaved bit sequence is generated by a bit interleaver before 
being modulated, e.g., Quadrature Amplitude Modulation (QAM), to a symbol sequence. Let iV r and N t 
denote the number of transmit and receive antennas respectively. Assume that S < mm{N t , N r } streams 
are transmitted at the same time. Hence, an S x 1 symbol vector is transmitted at the A;th time instant. 

The MIMO flat fading channel is assumed to be quasi-static and known by both the transmitter and 
the receiver, which is given by H e C NrXNt , where C stands for the set of complex numbers. Then, the 
beamforming matrices are determined by SVD of H, i.e., H = XJAV H , where U and V are unitary 
matrices, and A is a diagonal matrix whose sth diagonal element, \ s G M + , is a singular value of H in 
decreasing order with s = 1, . . . , S, where IR + denotes the set of positive real numbers. When S streams 
are transmitted at the same time, the first S columns of U and V, i.e., Us and V 5 , are chosen to be used 
as beamforming matrices at the receiver and transmitter subcarrier, respectively. Therefore, the system 
input-output relation at the A;th time instant is 

y kjS = \ s x k , s + n kjS , (2) 

with s — 1, . . . , S, where y ks and x kjS are the sth element of the S x 1 received symbol vector y fc and 
the transmitted symbol vector X& respectively, and n k , s is the additive white complex Gaussian noise with 
zero mean. 

The location of the coded bit c k / within the transmitted symbol is denoted as k' — > (k,s,j), which 
means that c k > is mapped onto the jth bit position on the label of x k)S . Let x denote the signal set of 
the modulation scheme, and let x{ denote a subset of \ whose labels have b E {0, 1} at the jth bit 
position. By using the location information and the input-output relation in ©, the receiver calculates the 
Maximum Likelihood (ML) bit metrics for c k i = b as 

A [y k>s , c fe /] = min \y k) „ - \ s x\ 2 . (3) 

Finally, the ML decoder, which applies the soft-input Viterbi decoding [9] to find a codeword c with the 
minimum sum weight and its corresponding information bit sequence b, uses the bit metrics calculated 
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by © and makes decisions according to the rule given by |[26ll as 



c = argminy>[y M ,c fc ,]. 

c. ' * 



c 

k' 



(4) 



With a properly designed bit interleaver, the maximum achievable diversity of BICMB is given by iflOll . 
ED as 



Based on ©, BICMB can achieve both full diversity and full multiplexing as long as the condition 
R r S < 1 is satisfied. 



If the channel is frequency selective, OFDM can be used to combat the ISI caused by multipath prop- 
agation [fT2l . OFDM is a Frequency Division Multiplexing (FDM) scheme, whose subcarrier frequencies 
are chosen to be orthogonal to each other in order to remove the Inter-Carrier Interference (ICI) among 
subcarriers and provide high spectral efficiency. 

The structure of OFDM is presented in Fig. [21 IFFT is used to multiplex the modulated symbols 
{x(l), . . . ,x(M)} at the transmitter, where x(m) denotes the modulated symbol transmitted by the mth 
subcarrier and M denotes the number of subcarriers. Hence, the output symbols of IFFT {5(1), . . . , x(M)} 
are calculated as 



where i = y/— 1. Then, Cyclic Prefix (CP) {x(M — L cp + 1), . . . , x(M)} is added before transmitting as 
guard interval, where L cp denotes the length of CP and L cp is no less than the number of channel taps 
L. The purpose of CP is to increase the OFDM symbol duration to eliminate ISI without breaking the 
orthogonality of transmitted symbols. 

At the receiver, CP is removed from the received symbols {y(M — L cp + 1), . . . , y(M),y(l), . . . , y(M)}. 
Then, FFT is used to derived the recovered symbols {y(l), . . . , y(M)} as 



D = (N r - \R C S] + l)(N t - \R C S] + 1). 



(5) 



B. OFDM 




(6) 



k=l 



M 




(7) 



k=l 
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Fig. 2. Structure of OFDM. 
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The frequency selective channel with L taps is given by h{l) E C with I = 1, . . . , L. Let 



h{m) = Y J W) 



27r(m— l)r^ 
g " MT 



1=1 



(8) 



denote the quasi-static flat fading channel observed at the mth subcarrier, where T denotes the sampling 
period and T\ indicates the /th tap delay ||27l . Then, the relation of x{m) and y(m) is given by 



y(m) = h(m)x(m) + n(m), 



(9) 



where n{m) is the additive white complex Gaussian noise with zero mean. As a result, symbol-by- symbol 
demodulation can be carried out for each subcarrier. 

OFDM is robust against ISI and has small sensitivity to time synchronization errors. It can be efficiently 
implemented using FFT and IFFT, and easily adapted to severe channel conditions without the compu- 
tational burden of time-domain equalization. It provides high spectral efficiency. Moreover, multipath 
diversity can be achieved by adding channel coding |fT3l , lfT4ll . The advantages of OFDM make it well- 
suited for broadband data transmission, and it has been selected as the air interface for all broadband 
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Fig. 3. Structure of BICMB-OFDM. 



wireless communication standards, i.e., WiFi, WiMax, and LTE [11511 . 



III. System Model of BICMB-OFDM 

BICMB-OFDM was proposed to achieve both spatial diversity and multipath diversity for MIMO 
frequency selective channels, ||20l , ETTl . 0. The structure of BICMB-OFDM is presented in Fig. [3j The 
channel coding, bit interleaver, and modulation of BICMB-OFDM are the same as BICMB introduced in 
section III-Al Assume that the modulated symbol sequence is transmitted through M subcarriers, and S 
streams are transmitted for each subcarrier at the same time. Hence, an S x 1 symbol vector x fc (m) is 
transmitted through the mth subcarrier at the kth time instant with m = 1, . . . , M. 

The MIMO frequency selective channel with L taps is assumed to be quasi-static and known by both 
the transmitter and the receiver, which is given by H(/) e QN r xN t w ^ / = 1, . . . , L. Let 

L 

. Inirn— 1)t^ 



H(m) = J2^-( l ) e ' i: 



(10) 



i=i 



denote the quasi-static flat fading MIMO channel observed at the mth subcarrier. Then, SVD beamforming 
is carried out for each subcarrier. The beamforming matrices at the mth subcarrier are determined by 
H(m) = U(m)A(m)V iT (m), where U(m) and V(m) are unitary matrices, and A(m) is the singular 
value matrix of H(m). When S streams are transmitted at the same time, the first S columns of U(m) 
and V(m), i.e., Ug(ra) and Vs(m), are chosen to be used as beamforming matrices at the receiver and 
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transmitter at the mth subcarrier, respectively. 

The multiplication with beamforming matrices is carried out at each subcarrier before executing the 
Inverse Fast Fourier Transform (IFFT) and adding CP at the transmitter, and after executing the Fast 
Fourier Transform (FFT) and removing CP at the receiver, respectively. The length of CP is L C p > L. 
Therefore, the system input-output relation for the mth subcarrier at the A;th time instant is 

yk,s( m ) = K(m)x ktS (m) + n k>s (m), (11) 

with s = 1, . . . , S, where y k ^ s {m) and x ktS (m) are the sth element of the 5x1 received symbol vector 
yfc(m) and the transmitted symbol vector x fc (m) respectively, and n k:S (m) is the additive white complex 
Gaussian noise with zero mean and variance iVo = N t /^ ll28l . with 7 denoting the received Signal-to- 
Noise Ratio (SNR) over all the receive antennas. Note that the total transmitted power is scaled by N t in 
order to make the received SNR 7. 

Similarly to BICMB introduced in section III-Al the location of the coded bit Cy within the transmitted 
symbol is denoted as k' — > (k,m,s,j), which means that cy is mapped onto the jth bit position on the 
label of x k)S (m). By using the location information and the input-output relation in (fTTI) . the receiver 
calculates the ML bit metrics for c k > = b G {0, 1} as 

A [yk, s (m),c k ,] = min \y k , s {m) - \ s (m)x\ 2 . (12) 
xex h< 

Finally, the ML decoder, which applies the soft-input Viterbi decoding to find a codeword c with the 
minimum sum weight and its corresponding information bit sequence b, uses the bit metrics calculated 
by (PT2l) and makes decisions as 

c = argmin A [y ktS (m), cy) . (13) 

c *— ' 

k> 

IV. Maximum Achievable Diversity of BICMB-OFDM 

The performance of BICMB-OFDM is dominated by the instantaneous Pairwise Error Probability (PEP) 
with the worst diversity. Based on the bit metrics in (fT2l) . the instantaneous PEP of BICMB-OFDM between 
the transmitted bit codeword c and the decoded bit codeword c is written as 



Pr (c — > c) = E [Pr (c — > c | H(m), Vra)] 



< E 



exp 



m) 



4N n 



(14) 



where dn denotes the Hamming distance between c and c, d min is the minimum Euclidean distance in 
the constellation, ^ fc , d ^ stands for the summation of the d H values corresponding to the different coded 
bits between the bit codewords. 

Define an M x S matrix A, whose element a ms denotes the number of distinct bits transmitting through 



the sth subchannel of the mth subcarrier for an error path, which implies £) m=1 J2s=i 



a. 



d H - Let aL 



denote the mth row of A. Note that the a-spectrum here is similar to BICMB in the case of flat fading 
channels introduced in IfTOll , [fTTTl . Then (PT41) is rewritten as 



Pr (c -> c) < E 
< E 



exp 
H exp 



47V 

4iVn 



(15) 



In this section, the maximum achievable diversity is derived by considering the simplest case where no 
correlation exists among subcarriers. Note that correlation among subcarriers has a negative effect on the 
performance, which will be discussed in Section |VJ 

Consider the case that different MIMO delay spread channels are uncorrelated and have equal power, 
and each element of each tap is statistically independent and modeled as a complex Gaussian random 
variable with zero mean and variance 1/L, i.e., 



E 



0, l?l>, 

i 1 = 1' 



(16) 



where h u>v (l) denotes the (u,v)th element in H(/). Further consider that the channel taps are separated 
by a constant sampling time, i.e., 77 = (I — 1)T in (flOl) . Then, the correlation between the mth subcarrier 
and the m'th subcarrier with m 7^ m' is given by 

E [h u>v (m)h* UjV (m')] 



,Je [KA m )K,v( m )] E [hu,v( m ')K,v( m ')] 





E 




Ef=iE 
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1 



1 - 



- 27r(m — m, , )L 

e~ l at 



(17) 



L 



1-e 



. 27r(m — m f ) 
* M 



Note that in the case of L = M, p = 0. This result implies that all subcarriers are uncorrelated. Note 
that h U) „ = [h u>v (0), . . . , h u>v (M)} is a complex normal random vector, where h u ^ v {m) denotes the (u, v)th 
element in H(m). This implies that the components of h UjV are independent if and only if they are 
uncorrelated [29 1. Therefore, based on (fTTT) . all subcarriers are independent in the case of L = M. In the 
following part of this section, this special case is considered. Although this special case is not practical, 
its diversity analysis provides the maximum achievable diversity for the practical case. The reason is that 
correlation among subcarriers for the practical case has a negative effect on performance, which will be 
discussed in Section |V] 

Since subcarriers are independent in this special case, the singular value matrices A(m) are independent. 
Therefore, (fT5l) is further rewritten as 



For each subcarrier, the terms inside the expectation in (fT8l can be upper bounded by employing the 
theorem proved in ll30l . which is given below. 

Theorem 1. Consider the largest S < min{N t , N r } eigenvalues /x s of the uncorrelated central N r x N t 
Wishart matrix / T371/ that are sorted in decreasing order, and a weight vector r\ = [rji, ■ ■ ■ , r]s] T with non- 
negative real elements. In the high SNR regime, an upper bound for the expression £'[exp(— 7 Y^ s =i Vs^s)], 
which is used in the diversity analysis of a number of MIMO systems, is 



where 7 is SNR, ( is a constant, r] min = min^.^o {Vi}i=i> an d $ is the index to the first non-zero element 
in the weight vector. 




in 



(18) 




(19) 



Proof: See El. 



□ 
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By applying Theorem \T\ to (fTST ), an upper bound of PEP is 



( 



) 



-D; 



Pr (c ->■ c) < ] [ C. 



, min OL rn,min 



■7 



(20) 



m,a m ^0 



with D m = (iV r — <5 m + 1) (iVt — <5 m + 1), where ot m ,min denotes the first non-zero element in a m , 5 m 
denotes the index of a m , m i n in a m , and ( m is a constant. Therefore, the diversity can be easily found from 
(1201 , which is 



Because the error path with the worst diversity order dominates the performance, the results of (1201 and 
(|2TI) show that the maximum achievable diversity of BICMB-OFDM is directly decided by the a-spectra. 
Note that the a-spectra are related with the bit interleaver and the trellis structure of the convolutional 
code, and they can be derived by a similar approach to BICMB in the case of flat fading channels presented 
in ifTOl . or by computer search. 

An example is provided here to show the relation between the a-spectra and the diversity. Consider the 
parameters N t = N r = S = L = M = 2. Assume that the R c = 1/2 convolutional code with generator 
polynomial (5, 7) in octal is employed, and the bit interleaver applies simple bit rotation, i.e., the sth bit 
of the mth S bits are transmitted through the sth subchannels at the mth subcarrier at one time instant. 
In this case, the dominant error path has the a-spectrum A = [01; 2 2], which implies that 8% — 2 and 
82 — 1. Hence, D\ = 1, D 2 = 4. Therefore, the maximum achievable diversity order is D = D% + D 2 = 5. 



A method to derive the a-spectra is illustrated by the following simple example. For this example, the 
system is composed of a 4-state 1/2 -rate convolutional encoder and a spatial de-multiplexer rotating with 
an order of a, b, c, and d which represent the four streams of transmission. Fig. [4] represents a trellis 
diagram of this convolutional encoder for one period at the steady state. Since a convolutional code is 
linear, the all-zeros codeword is assumed to be the input to the encoder. To find a transfer function of a 
convolutional code and a spatial de-multiplexer, the branches are labeled as a combination of a^, b^ b , c^ c , 
and d^ d , where the exponent denotes the number of usage of each subchannel which contributes to detecting 
the wrong branch by the detector. Additionally, Z^ z , whose exponent satisfies 0z = Pa + 0b + Pc + fid, 
is included to get the relationship between the Hamming distance djj and a-spectrum of an error event. 




D 



(21) 



A. a-spectra 
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ab X/X. cd 




'13 "23 

Fig. 4. Trellis of 4-state 1/2-rate convolutional code with 4 streams 
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Furthermore, the non-zero states are symbolically labeled from X n to X 23 as in Fig. HI while the zero 
state is labeled as Xi if branches split and X Q if branches merge, also as shown in Fig. HI 

Define x = [In, X 12 , X 13 , X 2 i, X 22 , X 23 \ . Then, the state equations are given by the matrix equation 



x = Fx + tXi 
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(22) 



We also get 



X = gx 



abZ 2 cdZ 2 



x. 



(23) 



The transfer function is represented in closed form by using the method in [|25l as 
T(a,b,c,d, Z) = g[I-F]- 1 t 

oo 

= gt + ^gF u t 



u=l 
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= Z 5 (a 2 b 2 d + bc 2 d 2 ) 

+ Z 6 (2a 2 bc 2 d + a 2 b 2 d 2 + b 2 c 2 d 2 ) 
+ Z 7 (a 2 b 3 c 2 + 2a 2 b 2 c 2 d + 2a 2 bc 2 d 2 + 

b\ 2 d 2 + aW + a 2 c 2 d?) 
+ Z\a% 2 c 2 + Aa 2 b 3 c 2 d + Aa 2 b 2 c 2 d 2 + 

b 4 c 2 d 2 + a 2 c 4 d 2 + Aa 2 bc 2 d 3 + a 2 b 2 d 4 ) + ■■■ (24) 

where [I — F] _1 can be expanded as I + F + F 2 + • ■ • through an infinite series of power of matrices. 

Consider BICMB-OFDM with parameters N t = N r = S = L = M = 2. Therefore, four streams are 
transmitted simultaneously. Assume that a and b are assigned to be the first and second streams transmitted 
by the first subcarrier respectively, and c and d are assigned to be the first and second streams transmitted 
by the second subcarrier respectively. The a-spectra can be figured out from the transfer function. For 
example, there are two error events with d h = 5 whose a-spectra are [2 2; 01] and [01; 2 2] respectively. 

This method can be applied to any /C-state k c /n c -rate convolutional code and BICMB-OFDM with M 
subcarriers and S streams transmitted at each subcarrier. If the spatial de-multiplexer is not a random 
switch for the whole packet, the period of the spatial de-multiplexer is an integer multiple of the Least 
Common Multiple (LCM) of n c and SM. Note that a period of the interleaver is restricted to correspond 
to an integer multiple of trellis sections. Define P = LCM(n c , SM) as the number of coded bits for a 
minimum period. Then, the dimension of the vector x is nP(fC — l)k c /n c where n is the integer multiple 
for a period of interest. 

B. Full Diversity Condition 

Based on the results of (|20l and (|2TT) . full diversity of N r N t L [[3] can be achieved if and only if all 
elements in the first column of the A matrix are non-zero, i.e., a m ,i ^ 0,Vm, for all error events. To 
meet such requirements of the a-spectra, the condition R C SL < 1 needs to be satisfied. In the following, 
the proof of the full diversity condition with the rate of the convolutional code is provided. 

Proof: To prove the necessity, assume that an information bit sequence b with length JR C SL is 
transmitted, then a bit sequence c m s containing J bits is transmitted at the sth subchannel of the mth 
subcarrier. If R C SL > 1, because the number of different codewords 2 JRcSL is larger than the number 
of different bit sequences c m s , 2 J , there always exists a pair of codewords result in the same c m>s . As a 
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result, the pairs of codewords with the same c m l cause full diversity loss. 

To prove the sufficiency, consider a bit interleaver employing simple rotation with the condition R C SL < 
1. In this case, all subchannels at each subcarrier could be assigned to one branch of the trellis structure 
of the convolutional code. Since the trellis of the convolutional can be designed such that the coded 
bits generated from the first branch splitting from the zero state are all error bits of an error event, each 
subchannel of all subcarriers could be used at least once for all error paths, which guarantees a m ,i 7^ 0, Vm. 
Therefore, full diversity can be achieved. 

This concludes the proof. □ 
The proof of the necessity above implies that in the case of R C SL > 1, there always exists an error 
path with no error bits transmitted through the first subchannel of a subcarrier. Therefore, full diversity 
cannot be achieved. In this case, the bit interleaver should be designed such that consecutive coded bits 
are transmitted over different subchannels of different subcarriers to provide the maximum achievable 
diversity, which depends on the a-spectra. 

V. Negative Effect of Subcarrier Correlation 

In practice, M is always much larger than L. In this case, correlation exists among subcarriers as shown 
in (fTTT ) and 021 . Hence, to calculate (fl"5l) , the joint Probability Density Function (PDF) of A(m)\ H (m) 
for all m satisfying a m ^ 0, which are eigenvalues of a set of correlated Wishart matrices OT1 . is 
required. However, this is an extremely difficult problem. The PDF of eigenvalues of two correlated 
Wishart matrices are given in 031 , 041 . 051 , which is already highly complicated. To the best of our 
knowledge, the joint PDF of eigenvalues of more than two correlated Wishart matrices is not available 
in the literature. The maximum diversity of an OFDM-MIMO system is known to be N r N t L $3$. In our 
case, however, a performance degradation caused by correlation is to be expected. Because, otherwise, 
the diversity can exceed the full diversity of N r N t L, which is a contradiction. In this section, the negative 
effect of correlation on the performance between two subcarriers is investigated to provide an intuitive 
insight. 

Consider an error path whose dn distinct bits between the bit codewords are all transmitted through 
two correlated subcarriers with correlation p, which could be the practical case. Define X = max(iV t , N r ) 
and Y = min(iV t , N r ). Let <fr = . . . , 0y] and <J> = [(pi, . . . , 0y] denote the ordered eigenvalues of 
the two correlated Wishart matrices HH H and HH fl , respectively. Note that 4> u = in (fl3T ). Let a = 
[ai, . . . , ay] and a = . . . , ay) denote the a-spectra of $ and $ respectively. Define p = \pi, . . . ,pw) 
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and p = [pi, . . . ,Pw) whose elements are the indices corresponding to non-zero elements in a and a, 
respectively, i.e, a Pw ^ and a p , ^ 0. Similarly, define q = [q u . . . , q Y -w] and q = . . . , q Y -w\ 
whose elements are the indices corresponding to zero elements in a and a, respectively, i.e, a Qw = and 



a d - =0. Define <fr D 



ypii • • • > Yp w i-> 
Therefore, the PEP in (fT5l) is written as 

Pr (c — > c) < E |^exp 

< E < exp 



'q\ i • • • i ~Tqu 



4AT 
w w 



. u;=i 



u)=l 



(25) 



with /i = {(^ min a m in) I (4iVo), where a min indicates the minimum element in a and a. To solve ((251) . the 
marginal PDF /(<& p , <3?p) is needed by calculating 



The joint PDF /($, $) is available in 03. OH, as 



/(*,*) d$ q d$ q . 



(26) 



/(#, <fr) = exp 
with the polynomial $) defined as 



1 y / 



(27) 



A(*,*) = tn(^-^)(^-^)] x K^4) (X " y)/2 /x-y(2Ve0^)l 



(28) 



M<1) 



where |-H" UiV | represents the determinant of the matrix with the (u, v)th element given by H u v , /„(•) is a 
modified Bessel function of order u which is given by 

2j+N 



In® = 



3=0 



j!(j + iV + l)! V2 



(29) 



and e ~ p 2 /(l — p 2 ) 2 - Because the exponent of p is related to the diversity, the constant appearing in the 
literature is ignored in (|28l ) for brevity. 

Since the eigenvalues of the Wishart matrices are positive and real, exp(j^(p u ) < 1 and exp(j^(f) u ) < 
1 in (T27T) . By applying Q u t e~ u &u < j^v t+1 and u l e~ u du = t\ to <fr q and $ q , the marginal PDF 
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/(3> p <I>p) in (|26l) is upper bounded as 



/(* P , * P ) < exp 



i-p 2 



X / 2 ($ p ,$ p 



p/' 



(30) 



where /2(3? p , ^ P ) is a polynomial corresponding to (l27l) . Then (|25l) is rewritten as 



Pr (c — > c) < 



o </o 



exp 



o </o 
ir 



ir 



/i + 



1 -p 2 



E^+E' 



, UI = 1 



jZ=1 



x/ 2 ($ pi $p)d$ p d$p. (31) 



Note that since /^(^p, ^p) is a polynomial, its multivariate terms can be integrated separately, and the 
term with the worst performance dominates the overall performance. To solve (|3TI) . a theorem provided 
in 11301 can be applied to integrate <3> p and <t> p independently for each multivariate term, which is given 
below. 

Theorem 2. For a multivariate term with variable (p w for w = 1, . . . ,W whose exponent, denoted by (3 W , 
is a non-negative integer, the multiple integration in a domain oo > <p\ > ■ ■ ■ > <pw > is 

„ , w 



OO P{pl 



JO 



w=l 



>W 



d0 1 = cM w+E »= i/ H 



(32) 



where ( is a constant. 
Proof: See [HI. 



□ 



Based on Theorem[2l the multivariate term with the smallest degree in /^(^p? 3*p) results in the smallest 
degree of (p + jz^)^ 1 , which dominates the overall performance. The smallest degree of /2(^ P , 3>p) is 
(X - pi + 1)(Y - pi + 1) + (X - pi + l)(y - p! + 1) - W - W, and the proof is provided in Appendix 
A. Therefore, (f3TT) is upper bounded by 

-D 



Pr (c -> c) < C 



d 2 ■ a 



1 



4X, 



-7 



I-P 2 



(33) 



with £> = (X - pi + 1)(Y - pi + 1) + (X - pi + 1)(Y - p~ x + 1), where C is a constant. 

The negative effect of subcarrier correlation is proved by (l33l) . When 7 — > 00, the diversity is the same 
as the uncorrelated case. However, on the practical range, the performance is degraded due to the term 
jzy, which is independent of SNR. Specifically, when the subcarrier correlation p is small, is also 



17 

relatively small, and its effect on the performance is not significant when the SNR is relatively large, and 
the uncorrelated case p = offers the performance upper bound. On the other hand, when p is large, 
is also relatively large compared to the SNR, then significant performance degradation could be caused, 
depending on the SNR. When p = 1, which means all the distinct bits of the error path are transmitted 
through one subcarrier, no multipath diversity is achieved, and the diversity equals BICMB in the case 
of flat fading channels introduced in ifTOl . ifTTI . which provides the performance lower bound. Note that 
the analysis in this section is not limited to equal power channel taps, and can also be applied to unequal 
power channel taps, non-constant sampling time, and other assumptions, which cause different subcarrier 
correlation. 

VI. Subcarrier Grouping 

The idea of subcarrier grouping technique is to transmit multiple streams of information through multiple 
group of subcarriers of OFDM. It was suggested for multi-user interference elimination [22J, Peak-to- 
Average Ratio (PAR) reduction E3l . and complexity reduction ll24ll . In this paper, the subcarrier correlation 
technique is applied to overcome the performance loss caused by subcarrier correlation. 

Note that p = in (fT71) when (m — m')L/M G Z, where Z denotes the set of integer numbers. 
This means that although correlation exists among subcarriers for L < M, some subcarriers could be 
uncorrelated if M/L e Z. In this case, there are M/L groups of L uncorrelated subcarriers. As a result, 
the subcarrier grouping technique is applied to transmit multiple streams of bit codewords through these 
M/L different groups of uncorrelated subcarriers, instead of transmitting one stream of the bit codeword 
through all the correlated subcarriers. As a result, the negative effect of subcarrier correlation is completely 
avoided, and the maximum achievable diversity is thereby achieved. For example, consider the case of 
L = 2 and M = 64. Then, the nth and the [u + 32)th subcarriers are uncorrelated for u — 1, . . . , 32. The 
subcarrier grouping technique can transmit 32 streams of bit codewords simultaneously through the 32 
groups of two uncorrelated subcarriers without performance degradation. Note that the best choice of the 
number of subcarriers in one group is L, since smaller choice results in less diversity while larger choice 
causes subcarrier correlation which degrades performance. 

Fig. [5] presents the structure of BICMB-OFDM with subcarrier grouping. In Fig. [5l Ti is a permutation 
matrix at the transmitter distributing modulated scalar symbols from different streams to their corre- 
sponding subcarriers, while T 2 = T^ 1 is a permutation matrix at the receiver distributing received scalar 
symbols of different subcarriers to their corresponding streams for decoding. Compared to BICMB-OFDM 
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Fig. 5. Structure of BICMB-OFDM with subcarrier grouping. 
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Fig. 6. Structure of BICMB-OFDM with subcarrier grouping in the frequency domain for one bit stream transmission. 



without subcarrier grouping, BICMB-OFDM with subcarrier grouping achieves better performance with 
the same transmission rate and decoding complexity. Note that the structure of BICMB-OFDM with 
subcarrier grouping can also be considered as Orthogonal Frequency-Division Multiple Access (OFDMA) 
|fT5l version of BICMB-OFDM. OFDMA is a multi-user version of the OFDM and it has been used in the 
mobility mode of WiMax ifTTll as well as the downlink of LTE [18J. The multiple access in OFDMA is 
achieved by assigning subsets of subcarriers to individual users, which is similar to the subcarrier grouping 
technique. As a result, with subcarrier grouping, BICMB-OFDM can provide multi-user compatibility. 

Fig. [6] presents the structure of BICMB-OFDM with subcarrier grouping in the frequency domain for one 
bit stream transmission. Note that Fig. [6] also presents the structure of BICMB-OFDM in the frequency 
domain when L = M. Therefore, the diversity analysis for L = M in Section [IV] can be applied to 
BICMB-OFDM with subcarrier grouping. As a result, the full diversity condition R C SL < 1 holds for 
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BICMB-OFDM with subcarrier grouping. In this paper, the number of employed subchannels by SVD 
for each subcarrier is assumed to be the same, which is S. However, they could be different in practice. 
In that case, the full diversity condition is R c J2t=i $u,i < 1 for u = 1, . . . , M/L, where S U) i denotes the 
number of employed subchannels by SVD for the Zth subcarrier of the wth group. 

Note that (TTTT) is derived under the assumption that all channel taps have the same power. When they 
have different power, there are no uncorrected subcarriers in general. However, some of them could 
have weak correlation. Therefore, the subcarrier grouping can still be applied to combat the diversity 
degradation, although it now can no longer fully recover the maximum achievable diversity because of 
subcarrier correlation. 

VII. Simulation Results 

To verify the diversity analysis, 2 x 2 M = 64 BICMB-OFDM with L = 2 and L = 4 using 4-QAM 
are considered for simulations. The number of employed subchannels for each subcarrier is assumed to 
be the same. The generator polynomials in octal for the convolutional codes with R c = 1/4 and R c = 1/2 
are (5,7,7,7), and (5,7) respectively, and the codes with R c = 2/3 and R c = 4/5 are punctured from 
the R c = 1/2 code ||25). The length of CP is L CP = 16. Each OFDM symbol has 4//s duration, of which 
0.8/is is CP. Equal and exponential power channel taps are considered. For the exponential channel model 
11361 , the ratios of non-negligible path power to the first path power are — 7dB, the mean excess delays 
are 30ns for L = 2 and 65ns for L = 4, respectively. The bit interleaver employs simple rotation. Note 
that simulations of 2 x 2 L = 2 and L = 4 BICMB-OFDM are shown in this section because the diversity 
values could be investigated explicitly through figures. 

Fig. shows the Bit Error Rate (BER) performance of 2 x 2 L = 2 M = 64 BICMB-OFDM 
employing subcarrier grouping over equal power channel taps with different S and R c . The A matrices 
that dominate the performance are provided in the figure. The diversity results of all curves equals the 
maximum achievable diversity orders derived from Section [TV] which is directly decided by the A matrices. 
Specifically, in the cases of S = 1, R c = 1/4 and R c = 2/3 codes, whose dominant A matrices are 
A = [2; 3] and A = [0; 5] respectively, achieve diversity values of 8 and 4 respectively. For S = 2, the 
codes with R c — 1/4, R c — 1/2, R c — 2/3, and R c — 4/5, whose dominant A matrices are A = [2 3; 3 3], 
A = [0 1; 2 2], A = [0 0; 2 2], and A = [0 0; 4] respectively, offer diversity of 8, 5, 2, and 1 respectively. 
Note that full diversity of 8 is achieved with the condition R C SL < 1. 
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Fig. 7. BER vs. SNR for 2x2L = 2M = 64 BICMB-OFDM with subcarrier grouping over equal power channel taps. 

Similarly, Fig. [8] shows the BER performance of2x2L = 4M = 64S' = l BICMB-OFDM employing 
subcarrier grouping over equal power channel taps with different R c . The A matrices that dominate the 
performance are provided in the figure. The diversity results of all curves equals the maximum achievable 
diversity orders derived from Section [IVj which is directly decided by the A matrices. Specifically, 
the codes with R c = 1/4, R c = 1/2, R c = 2/3, and R c = 4/5, whose dominant A matrices are 
A = [2; 3; 3; 3], A = [0; 1; 2; 2], A = [0; 0; 2; 2], and A = [0; 0; 0; 4] respectively, offer diversity of 
16, 12, 8, and 4 respectively. Note that full diversity of 16 is achieved with the condition R C SL < 1. 

Fig. [9] shows the BER performance of examined PEPs in (fl"5l) with S = 2, where the simplest case of 
an error event with du = 2 is examined for two subcarriers with different correlation coefficient p, which 
is derived from the 2x2L = 2M = 64 BICMB-OFDM over equal power channel taps. The figure shows 
that when p = 0, which implies the two subcarriers are uncorrelated, A = [1 0; 1 0] and A = [1 0; 1] 
offer diversity of 8 and 5 respectively. On the other hand, when p ^ 0, performance degradation is caused 
by subcarrier correlation, and stronger correlation results in worse performance loss. When p = 1, which 
means the d H = 2 distinct bits are transmitted through only one subcarrier and no multipath diversity is 
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Fig. 8. BER vs. SNR for2x2L = 4M = 64S' = l BICMB-OFDM with subcarrier grouping over equal power channel taps. 

achieved, both A = [2 0; 0] and A = [1 1; 0] provide diversity of 4. The results are consistent with the 
analysis provided in Section |Vl and they show the negative effect of subcarrier correlation on performance. 

Fig. M shows the BER performance of 2 x 2 I = 2 M = 64 5 = 1 fi c = 1/2 BICMB-OFDM with 
and without subcarrier grouping over equal and exponential power channel taps. In the figure, w/ and 
w/o denote with and without respectively, while SG denotes subcarrier grouping. The results show that 
the subcarrier grouping technique can combat the performance loss caused by subcarrier correlation for 
both equal and exponential power channel taps. As discussed in Section [VH the maximum achievable 
diversity of 8 is provided by employing subcarrier grouping for equal power channel taps, since there is 
no subcarrier correlation. As for the case of exponential power channel taps, subcarrier grouping cannot 
fully recover the performance loss because subcarrier correlation still exists. 

Fig. ED] shows the correlation p of two subcarriers with different separation for 2x2L = 2M = 64 
BICMB-OFDM over equal and exponential power channel taps. The figure shows that the channel with 
exponential power taps causes stronger subcarrier correlation than equal power taps, which results in 
worse performance as shown in Fig. [TOj 
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Fig. 9. BER vs. SNR for examined PEPs of two subcarriers with different correlation for 2 X 2 L = 2 M = 64 S = 2 BICMB-OFDM 
over equal power channel taps. 



Similarly to Fig. EE Fig. [12] shows the BER performance of 2 x 2 L = 4 M = 64 S = I R c = 1/4 
BICMB-OFDM with and without subcarrier grouping over equal and exponential power channel taps. 
The results show that the negative effect of subcarrier correlation on performance can be overcome by the 
subcarrier grouping technique for both equal and exponential power channel taps. For equal power channel 
taps, the maximum achievable diversity of 16 is achieved by applying subcarrier grouping as discussed in 
Section[VT]because subcarrier correlation is totally removed. On the other hand, since subcarrier correlation 
is not totally reduced in the case of exponential power channel taps, subcarrier grouping cannot fully restore 
the performance degradation. 

Similarly to Fig. [TTJ Fig. [13] shows the correlation p of two subcarriers with different separation for 
2x2L = 4M = 64 BICMB-OFDM over equal and exponential power channel taps. The figure shows 
that the subcarrier correlation of exponential power channel taps is larger than equal power channel taps 
so that it achieves worse performance as shown in Fig. f 
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Fig. 10. BER vs. SNR for 2 x 2 L = 2 M = 64 S = 1 iic = 1/2 BICMB-OFDM with and without subcarrier grouping over equal and 
exponential power channel taps. 



A. Discussions 

Fig. |9] verifies the negative effect of subcarrier correlation as analyzed in Section |Vj while Fig. [TOland 
Fig. [[2] verify the advantage of BICMB-OFDM with subcarrier grouping over BICMB-OFDM without 
subcarrier grouping. Other than the subcarrier grouping technique, focus can be drawn on the design of 
the bit interleaver to combat the subcarrier correlation. As illustrated in Section IIV-AI the bit interleaver 
is directly related to the a-spectra which reflect the subcarrier distribution of the distinct bits for each 
error path and thereby determine the diversity. In fact, the performance degradation results from the error 
paths whose distinct bits are transmitted by correlated subcarriers. Therefore, negative effect of subcarrier 
correlation can be reduced by a properly designed bit interleaver. However, such an interleaver only finds 
the better arrangement of subcarriers to lighten the negative effect of subcarrier correlation without actual 
reduction achieved by subcarrier grouping. As a result, the subcarrier grouping technique is apparently a 
better choice because it does not only achieve better performance but also provides multi-user compatibility 
as explained in Section |VIl Therefore, the design of the bit interleaver to combat subcarrier correlation 
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Fig. 11. Correlation vs. subcarrier separation for 2x2L = 2M = 64 BICMB-OFDM over equal and exponential power channel taps. 



is not considered in this paper. 

In this paper, FFT-based OFDM is considered in which CP is added to combat ISI. However, adding 
CP results in data rate loss. In order to avoid the extra CP data that reduce bandwidth, Discrete Wavelet 
Transform (DWT) and Inverse DWT (IDWT) has been considered as alternative platforms for OFDM 
to replace IFFT and FFT 03, (31, ED. DWT-based OFDM is overall the same as FFT-based OFDM, 
and the only difference is in the OFDM modulator and demodulator. The DWT-based OFDM employs 
Low Pass Filter (LPF) and High Pass Filter (HPF) operating as quadrature mirror filters satisfying perfect 
reconstruction and orthonormal bases properties. The transform uses filter coefficients as approximate and 
detail in LPF and HPF respectively. The main advantage of DWT-based OFDM is that it does not need 
adding CP to combat ISI, and thereby provides higher data rate than FFT-based OFDM. However, the 
DWT-based OFDM requires more complicated equalization ||40l , [14111 . [|42l . than the simple frequency- 
domain single-tap equalization of FFT-based OFDM for frequency selective MIMO channels. As a result, 
the corresponding quasi-static flat fading MIMO channel for each frequency subcarrier in (flOl) may not 
be valid so that the SVD beamforming cannot be employed. Therefore, the FFT-based BICMB-OFDM 
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Fig. 12. BER vs. SNR for 2 x 2 L = 4 M = 64 S = 1 i? c = 1/4 BICMB-OFDM with and without subcarrier grouping over equal and 
exponential power channel taps. 



proposed in this paper may not be applied to DWT-based OFDM. Nevertheless, if a single-tap equalization 
similar to © becomes available for DWT-based OFDM, then the BICMB-OFDM technique and its 
diversity analysis proposed in this paper become valid. The only difference is the subcarrier correlation so 
that a different selection needs to be considered for subcarrier grouping. The discussion on the existence 
of single-tap equalization technique for DWT-based OFDM is beyond the scope of this paper, so it is left 
for future considerations. 

Fig. Hand Fig. [8] verify the relation between the diversity and the a- spectra as well as the full diversity 
condition R C SL < 1 derived in Section [IV] for BICMB-OFDM with subcarrier grouping. The full diversity 
condition implies that if the number of streams S transmitted at each subcarrier increases, the code rate R c 
may have to decrease in order to keep full diversity. As a result, increasing the number of parallel streams 
may not fully improve the total transmission rate, which is a similar issue to the full diversity condition 
R C S < 1 of BICMB for flat fading MIMO channels introduced in OH, 0J]. In fact, for flat fading MIMO 
channels, other than channel coding, the constellation precoding technique has been incorporated with both 
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Fig. 13. Correlation vs. subcarrier separation for 2 x 2 L = 4 M — 64 BICMB-OFDM over equal and exponential power channel taps. 

uncoded and coded SVD beamforming to achieve full diversity and full multiplexing simultaneously, with 
the trade-off of a higher decoding complexity [|43l . PBl . [|45l . Il46ll . Specifically, in the uncoded case, 
full diversity requires that all streams are precoded. On the other hand, for the coded case, which is 
BICMB, even without the condition R C S < 1, other than full precoding, partial precoding with lower 
decoding complexity than full precoding could also achieve both full diversity and full multiplexing with 
the properly designed combination of the convolutional code, the bit interleaver, and the constellation 
precoder. Moreover, Perfect Space-Time Block Codes (PSTBCs) ll47l . which have the properties of full 
rate, full diversity, uniform average transmitted energy per antenna, good shaping of the constellation, and 
nonvanishing constant minimum determinant for increasing spectral efficiency which offers high coding 
gain, have been considered as an alternative scheme to replace the constellation precoding technique for 
both uncoded and coded SVD beamforming with constellation precoding in order to reduce the decoding 
complexity in dimensions 2 and 4 while achieve almost the same performance [|48l , 11491 , ||50l . Since 
these techniques have been successfully solved the restricted full diversity condition issue of R C S < 1 
for BICMB in the case of flat fading MIMO channels, it may be possible to incorporate these techniques 
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with BICMB-OFDM so that its full diversity condition is not restricted to R C SL < 1 for frequency 
selective MIMO channels. However, the design criteria and diversity analysis cannot be generalized in a 
straightforward manner because of the increased system complexity so that they are considered as future 
work. 

VIII. Conclusions 

BICMB-OFDM combines MIMO and OFDM to achieve both spatial diversity and multipath diversity 
for frequency selective MIMO channels so that it can be an important technique for broadband wireless 
communication. In this paper, the diversity analysis of BICMB-OFDM is carried out. As a result, the 
maximum achievable diversity is derived and the important a-spectra directly determining the diversity 
is introduced, providing important insights of BICMB-OFDM. According to the analysis, a sufficient and 
necessary condition R C SL < 1 for achieving full diversity is proved, which is very important for practical 
design. In addition, the negative effect of subcarrier correlation on the performance in the practical case is 
investigated, and subcarrier grouping is employed to overcome the performance degradation and provide 
multi-user compatibility. Furthermore, it may be possible to employ precoding techniques so that the full 
diversity condition of BICMB-OFDM is not restricted to R C SL < 1, and its diversity analysis and design 
criteria are considered as future work. 

Appendix A 
Proof of the Smallest Degree of / 2 (<& p , <J> p ) 

The polynomial fzi^p, 3>p) in (|30l) corresponds to dTTT) . Because i^e~ u du < ^-f t+1 and J °° u l e~ u du = 
t\, the smallest degree of /^(^p, ^p) is related to the polynomial <&) in (T27T) . which is given by 

(|28T) and can be rewritten as 

Y Y 

*) = e {x - Y ^ 2 [H(<p u - <j> v ){4> u - ^)][J](0 M ^) x - y ]|/ X -y(e</) u 1 ,)|, (34) 

u<v u=l 

where 

Note that only the multivariate term of /i($, $) determining the smallest degree of /^(^p, ^p) needs 
to be considered. The dominant term of fi(<f>,<$>) is the one with the smallest degree and the largest 
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eigenvalues, which depends on the dominant term of n«<i;(0u — 4>v){4>u — 4>v) an d the dominant term of 
\Ix-Y{e<t>u4>v)\- The dominant term of Yll <v ((f> u ~ <Pv)(<Pu - 4> v ) is \~^ a= i{<Pu<Pu) Y ~ u ■ The dominant term 
of \Ix-y(e(pu4>v)\ is C 11^=1 (^«^«) y u wnere C is a constant, and the proof is provided in Appendix B. 
Therefore, the dominant term in /i($, ignoring the constant, is given by 

Y 
u=l 

The degree of /i($, is 

6 h =2[Y(Y-l) + Y(X-u)}. (37) 

After integration of ([26]), the factor <Pu +Y ~ 2u and me factor Wu=i 4>u +Y ~ 2u of vanish 

because f °° u t e~ u du = t\. Therefore, 

6 V anished = (pi ~ l)(X + Y - Pl ) + (p x - 1)(X + Y - px). (38) 

Meanwhile, the eigenvalues with q u > px and (pg u with g u > px result in increased degree because 
Jq u t e~ u du < ^j-t> <+1 . Therefore, 

^dded = 2Y - W - W - V x - Pi + 2. (39) 

As a result, the smallest degree of /^("^p, *£ P ) is 

& ^ $ vanished ~\~ $ added 

= (X - p! + 1)(Y - px + 1) + (X - px + l)(y - p! + 1) - W - W. (40) 

Appendix B 

PROOF OF THE DOMINANT TERM OF \Ix-Y{e<j>u<f>v)\ 

When Y = 1, 

l^'I^S^i (41) 

and the dominant term is 1/X\. 
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When Y — 2, 



\I X -y{^u^v)\ = Ix-Y(e4>i4>i)Ix-Y(e4>2<t>2) - ^x-y(e0i02)^x-r(e020O 



E 



1)! ^ ? !(?+X-l)! 



, ./!«./ • -V • 1)! ^ ./!(./ • .V- 1)! ^ .,!(./ • .Y 



1)! ? 7! 



2 2 



u=i u=i 



EE 

.3=0 k>j 



j\(J + X - 1)\ k\(k + X - 1)\ 



and the dominant term is e0i0i/ [A"!(X — 1)!]. 
When Y > 3, 

y y 



^201; 



3=0 



jl(j + X-l)\ 



(42) 



|/x-y(e^«( 



11=1 v=l 



' Y oo 

nE 

.k=ljk=0 



\:ik 



j k \( Jk + X-Y + l)\ 



(43) 



3k<3k+l 

where u k = [(u + k — 2) mod Y) + 1 and v k = [(v + k — 2) mod Y) + 1, and the dominant term is 
CUl =1 (<PM Y ~ k with C = Il Y k=1 £ Y - k / l(Y - k)\(X - k + 1)!]. 
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